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Dependence of a/r c (inverse Sommer parameter in units of lattice spacing a) on am q (quark mass in lattice 
unit) has been observed in all lattice QCD simulations with sea quarks including the ones with improved actions. 
How much of this dependence is a scaling violation has remained an intriguing question. Our approach has been 
to investigate the issue with an action with known lattice artifacts, i.e., the standard Wilson quark and gauge 
action with j8 = 5.6 and 2 degenerate flavors of sea quarks on 16 3 x 32 lattices. In order to study in detail the 
sea quark mass dependence, measurements are carried out at eight values of the Wilson hopping parameter K 
in the range 0.156 - 0.158 corresponding to PCAC quark mass values am q from about 0.07 to below 0.015. We 
analyze the static potential by fitting to the familiar phenomenological form and extract a/r c . Though scaling 
violations may indeed be present for relatively large am q , a consistent scenario at sufficiently small am q seems 
to emerge in the mass-independent scheme where for a fixed j3, 1 /tq and -/a have linear dependence on m q as 
physical effects similar to the quark mass dependence of the rho mass. We present evidence for this scenario 
and accordingly extract the lattice scale a by chiral extrapolation to the physical point. Care has been exercised 
to determine optimal values of all fitting parameters and accuracy of the chiral extrapolation. An independent 
determination of the scale a by chiral extrapolation of the rho mass is consistent with the scale obtained above 
(a = 0.08041(12)(77) fm, a~ l = 2.454(4)(23) GeV). 

PACS numbers: 02.70.Uu, ll.10.Gh, ll.10.Kk, 11.15.Ha 



I. INTRODUCTION 



An accurate determination of the lattice scale is mandatory for comparing lattice observables with their continuum counter- 
parts. While the determination of lattice scale is conceptually simple in the quenched approximation of lattice QCD, simulations 
with dynamical quarks bring forth some unavoidable complications. It is of interest to study whether some of the complications 
are due to lattice artifacts or they throw light on physical aspects of lattice QCD. 

In the early days of lattice QCD, string tension provided a method to set the scale. If the color electric flux emanating from 
the quarks are squeezed into one dimension (a string-like configuaration), such a gauge field configuartion will have constant 
energy per unit length i.e., E = or. The energy density a is called the string tension. For the quenched theory, string tension 
is a fundamental quantity. For a review, see, Ref. |Tj. In lattice gauge theory, string tension is defined as the asymptotic value 
of the force F(r) between a pair of static quark and anti-quark [2] at separation r, i.e., string tension a = limit, ■^ a> F(r), The 
limiting value may not be easy to extract from the lattice data since statistical errors on the force (which is extracted from the 
large Euclidean time exponential fall-off of the expectation value of the Wilson loop) increases with the distance. 

In the presence of dynamical quarks, the string may break at large r and the definition of string tension as the asymptotic 
value of the force between static sources may not be useful. However, if we define the string tension to be the coefficient of the 
linearly rising potential (which in principle exists at all quark antiquark separations within the hadron i.e., for less than a fermi 
(fm)), then this notion can survive the presence of dynamical quarks and may be used for practical purposes. If we use the string 
tension to set the scale, we need to know the value of the string tension in physical dimensions and there exists a fair degree of 
uncertainty about this value J2] [3] ffl [5] |6] |7) . In addition, it was also noticed with the introduction of sea quarks that the string 
tension in lattice units has a dependence on the quark mass in lattice units. 

In Ref. 1 8 1 Sommer introduced the method of determining the lattice spacing through distance scales r c derived from the poten- 
tial between a static quark-antiquark pair, e.g. , using phenomenologically reasonably well-known information of rQ,F(ro)=1.65 
where ro = 0.49 fm. Also see Ref. [9|. The Sommer scale ro was originally conceived as a bosonic observable and was expected 
to be independent of the sea quarks. However, all lattice QCD simulations with dynamical quarks employing a variety of actions 
(including improved quark and gauge actions) have shown that the ratio ro/a, i.e., the Sommer scale in units of the lattice spacing 
a depends significantly on am q , the quark mass in lattice units. 
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In this work we carry out a detailed and careful investigation of the static potential and the ratio a/r c , where r c is a Sommer- 
type scale, in a lattice QCD simulation with standard Wilson gauge and quark actions on 16 3 x 32 lattices at a single gauge 
coupling parameter j3 = 6/g 2 = 5.6 with a large set (eight values) of the fermionic hopping parameter K for the sea quarks. 
Use of this large set of sea quark masses has helped us to identify, in terms of a certain parameterization of the static potential, 
the scale-violating part of the dependence of a/r c on am q . Our numerical data support the interpretation that for small enough 
am q (< 0.035) the dependence of a/r c on am q is a physical effect. 

In a mass-independent scheme, something that follows quite naturally as explained later in this paper, the scale a is then 
obtained by chirally extrapolating a/r c to the physical point. Since no other information is available on the dependence of a/r c 
on am q other than from numerical simulations, this extrapolation needs to be done with care to exclude large uncertainties from 

the extrapolation. Ignoring possible uncertainties of the values of rf 1 (r c at the physical point), we find accurate values of the 
scale determined this way with about 1 % error. 

Independent of the determination from the static potential and the Sommer scale, we have also determined the scale a from 
chiral extrapolation of am p , the rho meson mass in lattice unit. The scale determined this way is consistent with the scale 
determined from the chiral extrapolation of a/r c , although with somewhat larger error bars (~ 2 — 2.5%). 

In our determination of the static potential and the subsequent analysis to obtain the ratio a/r c , we have exercised utmost care 
in determining the fit range to determine the static potential, the fit range to determine the parameters of the static potential and 
also the optimum smearing levels to be used for the gauge configurations. In addition, for the chiral extrapolations of ao 1 / 2 and 
a/r c , we have first used am q and then (am K ) 2 as the chiral regulator to double check the reliability of the extrapolation. We have 
preferred (am K ) 2 to (r c m n ) 2 for chiral extrapolation to the physical point because r c itself has a chiral dependence. We have used 
all possible cross-checking of different determinations for consistency of our results and we present evidence for consistency in 
this paper. 

Our results are based on accurate determinations of the parameters of the static potential and reliable chiral extrapolations. 
We believe that our qualitative conclusions, if not also the quantitative conclusions, are independent of the numerical details like 
the fit-ranges, smearing levels etc. 

In an earlier paper [ 10 1, at the same set of parameters with the same action we have determined, using gaussian smearing both 
on source and sink, the pion and the rho masses, their decay constants and the PCAC quark mass. We have used most of those 
results in this paper. For details, please see Ref. ifTOl . 

Before we end this section, we would like to point out that there exists in the literature a variety of other methods to determine 
the lattice scale ifTTI [Til [El [T4l [El . 

This paper is organized as follows. Section [II] contains a summary of the simulation. In Sec. Ill measurement of the Wilson 
loop and the extraction of the static potential are given and Sec. |IV] describes the fit of the static potential. Sommer scale is 
described in Sec. IV] Dependence of various parameters of the potential on the lattice quark mass am q is detailed in Sec. VI and 
the interpretation or the am q dependence is presented in Sec. VII Section VIII presents the extrapolation of the data to chiral 
and physical points and the determination of the lattice scale at the physical point. Se c. |IX| presents estimates of the physical 
string tension. Implications from weak coupling perturbation theory are explored in Sec. |X| Finally, Sec. [Xljdiscusses the salient 
features and the conclusions. 



H. SIMULATION 

We have used unimproved Wilson gauge and fermion actions at a single gauge coupling given by j3 = 6/g 2 = 5.6 and two 
flavors of degenerate sea quarks on 16 3 x 32 lattices. Our choice of the gauge coupling is motivated by the requirements of a 
reasonably small lattice spacing so that the results of the simulation are not significantly contaminated by scaling vilolations. 
As mentioned already in the Introduction, lattice QCD simulations in the past using various gauge and fermion actions have 
observed significant sea quark mass dependence on off 1 / 2 and a/r c . In order to study this sea quark mass dependence in detail, 
we have generated gauge configurations and performed measurements at a large set of sea quark masses corresponding to eight 
values of the Wilson hopping parameter k = 0.156, 0.1565, 0.15675, 0.157, 0.15725, 0.1575, 0.15775 and 0.158. 

At each K we have generated 5000 equilibrated configurations with the standard HMC algorithm (with even-odd pre- 
conditioned Conjugate Gradient for inversion of M^M, M being the fermion matrix) and performed Wilson loop measurements 
separated by 25 configurations. Details on the simulation and autocorrelation times can be found in Ref. iflOl . 

In Ref. [ 10 1 we have presented a detailed study of the pion and the rho mass, their decay constants and PCAC quark mass using 
the same set of simulation parameters. We used gaussian smearing on both source and sink and investigated systematic effects 
on the pion mass and the decay constant using different types of correlators (PR AA, AP and PA, where P and A respectively are 
pseudoscalar and axial vector densities). PCAC quark masses were determined using different pion correlators. We observed 
finite size (FS) effects on our lowest pion masses. However, interestingly the effect was different for different operators, e.g., 
pion mass from the AA correlator at K = 0.15775 had negligible FS effect while the same computed from the PP correlator had 
significant FS effect. At K = 0.158 (our largest value of fc), pion masses computed from either operator had significant FS effect, 
but the effect was less in the AA correlator. 
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In this paper, we shall study the am q dependence of a/r c and other quantities. Eventually the pion masses in lattice units are 
used for chiral extrapolation of quantities derived from the static potential to the physical point. We use all results of quark and 
pion mass obtained in Ref. ifTUl in this paper. We use similar notation and convention as developed in ifTOl , e.g., lattice quark 
masses are denoted arn^ or am q p depending on the correlator used and each of these quark masses further depend on whether 
the pion mass used in the determination of the quark mass was taken from the PP or the AA correlator. 

Unless otherwise stated, all errors in this paper presented in data (in text or tables) or shown in figures are single-omission 
jackknife errors computed from 200 jackknife bins. 



III. EXTRACTION OF THE STATIC POTENTIAL FROM WILSON LOOPS 

For sufficiently large times T, the asymptotic behavior of the expectation value of the Wilson loop W(R,T) (R being the 
spatial separation) is given by (W(R,T)) = C(R) exp[— V(R)T~\ , where V(R) is the potential between a pair of static quark and 
antiquark. The coefficient C(R) is the ground state overlap. In order to determine the static potential reliably, it is important to 
have a large ground state overlap in the measurement of the Wilson loop. 

The Wilson loop is a gauge-invariant quantity and as such gauge-fixing the gauge configurations is not necessary for its 
measurement, although we do it anyway. After gauge fixing to temporal gauge [ 16 1, APE smearing 1 3 1 is performed on the link 
fields. Smearing gets rid of short distance fluctuations, helps in reducing higher states' contamination and increase overlap with 
physical states of interest. For a very early discussion of the need for smearing, see Ref. ifTTl . Smearing of gauge fields in a 
fermion action is also expected to have several other advantages lfl"8ll . 

APE smearing is performed as follows: 

U i {x)^U' i {x) = {l-c)U i {x)+ C - £ Ui(x), wheTeUi{x)=Uj(x + i)UJ{x+j)u](x), (3.1) 

staples 

followed by projection back to 51/(3). The parameter c is the relative strength of the smearing and we have chosen c = 4/(e + 4) 
with £=2.5. 

We measured the Wilson loops (W(R, T)) with temporal extents up to T = 16 and spatial separations up to R = y3 x 8. 
A reasonable estimate of the static potential aV (R) is obtained by the plateau reached at large T of the effective potential 

aV eS (R, T) = In , ( W ( RlT ^ . (3.2) 

Fig. [T]shows, at K = 0.15775 and smearing level 30, aV e ff as a function of T for a host of values of R ranging from R = 1 to 
R = 7. Generally plateaux are observed in this figure starting from R = 3. As T grows for a given R, or as R grows, the data get 
noisy. 

At each K the optimum level of smearing (which in the case of K = 0.15775 is 30 for the data shown in Fig. [TJ is obtained 
by comparing the ground state overlap C(R) for different smearing levels. This is what is shown in Fig. [2] At large values of R, 
the ground state overlap increases as the smearing level increases (upto a certain smearing level), but at small R (as shown in the 
inset) the reverse is generally true. In order to have the optimum ground state overlap we have chosen 30 as the smearing level 
at this K. For accurate analysis of aV (R) in a given range of R, in terms of parameters sensitive to either small R or large R, it is 
important to choose the smearing level which gives optimum ground state overlap throughout the range of R used in the analysis 
of aV(R). 

The optimum smearing levels used in this paper are 25 (for K = 0.156, 0.1565, 0.15675, 0.157, 0.15725) and 30 (for 
K = 0.1575,0.15775,0.158). 

For each value of R, we determine V(R) by a single exponential fit in the T range [r m j„, r max ] = [3,4], [3,5] and [4,5]. The 
single exponential fitting ansatz and the fitting ranges are well justified by the plateaux in Fig. [T] Eventually we have chosen the 
range [3,4] for the final analysis, as discussed later. 



IV. FIT OF THE STATIC POTENTIAL 



Phenomenologically lfT9ll . the potential V between a static quark and antiquark is parameterized as follows: V(r) = Vq + 
a r + j where a is the string tension which has the dimension of mass 2 . In lattice units, we have aV(r) = oVq + a 2 a ~ + a". 
Writing r = Ra and a = a/a 2 , we get aV(R) = aVo + &R+ f . 

After incorporating the correction for the finite lattice using the expression for the perturbative lattice Coulomb potential 
MM 
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(4.1) 




FIG. 1: The effective static potential aV e ff(R, T) as a function of T for different values of R at smear level 30 and K=0. 15775. 
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FIG. 2: Ground state overlap C(R) versus R for K=0. 15775 as a function of the smearing level (smJev). The inset shows the enlarged view of 
the small R region. The dashed line connecting the data points at sm_lev = 30 shows the optimum nature of this choice both at large and small 
R region. 



the parameterization of the potential on the lattice reads 

a /hi A 

" T, (4-2) 



aV(R) = aV + OR - - - <5rot 
K 



where 5rot is the coefficient of the correction term. The measured static potential is fit to the formula in Eq. (4.2 1 which 
corrects the lattice data for the lattice artifacts in the Coulomb potential. The first three terms of Eq. ( |4.2[ ) now gives the 
continuum potential (i.e., without lattice artifacts). 

There are a few comments on the 1 /R terms proportional to a and 5rot: (i) the lattice version of 1 /R emerges out of fourier 
transforming the gluon propagator 1/q 2 in a finite box for 1-gluon exchange interaction between a pair of static quark and 
antiquark and as such the parameter a has the interpretation of being proportional to the strong coupling constant, (ii) the 
difference ([1/1?] — 1//?) between the lattice version and the continuum version is never negligible on a finite lattice. As a result 
the correction is never very small as evident in Fig. [3] (iii) the parameter a is expected to run with R at these intermediate length 
scales, (iv) we can then only estimate an average a over the values of R where the static potetial is fit, (v) perturbative running 
is generally applicable at scales > 2 GeV which translates into R < 1 in our case. 

Fig. [3] shows in 2 plots at K = 0.1575 the uncorrected potential as obtained from our numerical simulation and the corrected 
(continuum) potential obtained by subtracting out the correction (5rot) term. The open and solid symbols are respectively the 
uncorrected and the corrected potential. The dotted straight lines are just joining the uncorrected points. The solid lines going 
through the corrected data represent the fit and by definition it should go through the corrected points within the /Grange of the 
fit which in this case is from R m [n = a/2 ~ 1.4 to R ma _ x = 3^/5 ~ 6.7. The purpose of the upper plot is to show that the kinks at 
small R region due to breaking of rotational invariance in our finite lattice disappear after the correction. The lower plot shows 
that even much beyond the fit-range the fit to the corrected potential data points are very good; although this is generally true at 
all K, in this case we have picked a particularly good example (i.e., at K = 0.1575). 

Because of noisy data at large R and also the exponential fall-off of the expectation value of the Wilson loop, the potential is 



poorly determined at large distances. As already noted above, the expression in Eq. (4.2 1 ignores the running of the coupling a 




FIG. 3: Potental with and without the correction for finite lattice for K= 0.1575. The upper panel shows the fitted range while the lower panel 
further shows the quality of the fit in regions of R beyond the fitted range. 
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FIG. 4: The parameter a for ff =0.1575 as a function of R max for different choices of R m i n for three different data sets corresponding to three 
T ranges. 



|. Hence it is advisable to fit the potential in as limited a range as possible [23 1. On the other hand, determination of a which 
has the interpretation of the string tension at large R becomes uncertain if R is not taken large enough. 

We fit the potential according to Eq. (4.2 1 in the range {R m in 5 Rmax} with R m { n = y/2, V3,2, V5, 2\/2 and R max = 
4^2,6,3^5, V3,7. 

Fig. |4] shows, at K = 0.1575 with APE smearing level 30, the fit value of the parameter a for different values of R mm and 
.Rmax- The figure also shows 3 differnt data sets corresponding to the 3 ranges [7m m ,7max] = [3, 4] , [3, 5] , [4, 5] used for the 
evaluation of the potential by a single exponential fit. The figure clearly shows that the errors for the fit with [r m ; n , T maK ] — [3,4] 
and R m [ n = \[2 are the minimum. In addition, for this choice, the data is also relatively independent of R max . For these reasons 
we have chosen the fit range for the V (R) evaluation as [7m m , T max ] — [3,4] and the fit range for fitting the potential to Eq. (4.2 1 

to be [/? min , ^max] = V2, 3\/5 . 

The parameter a is always determined better at smaller R m i n . We also notice in Fig. [4] that the value of a grows larger as 
^min increases, a fact consistent with the expected running of a. The results of the fit, i.e., the values of the parameters a, 6 
etc. are never fully invariant under the change of R m in- As already pointed above, a is sensitive to /? m in an d as a result the other 
parameters of the correlated fit follow suit. The parameter a is only approximately independent of R mm for Rmm > 2, but as far 
as any reliable determination of a is concerned, the region R m \ B > 2 is not reliable. 

Fig. [5] which plots a as a function of R m [ n for the three different fitting range [7min;7max] clearly shows that a depends on 
^min an d the uncertainty in a increases with increasing R m in- Based on this figure alone, perhaps the range [7mi n , 7max] = [3,5] 
is equally good as the range [3,4], but overall for all data the errors are smaller for our chosen range [3,4]. 

At this point, let us also point out that a depends, albeit quite weakly, also on the smearing level as more and more smearing 
progressively cuts out high frequency modes. We have already explained, based on Fig. |2j how we have determined the level of 
smearing at each K. 

The results of the fit to Eq. (4.2 1, as explained in detail above, depend to some extent on various parameters related to the 
fitting procedure (r m i n , r max , Rmin, Rmax) a nd a lso on the smearing level. We have found that our final conclusions regarding the 
Sommer parameter, its dependence on the quark mass and the lattice scale are not very sensitive to change of these parameters 
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FIG. 5: The parameter a for K =0.1575 as a function of R m \ n for three different T ranges. For each T range, the behavior for five choices of 
R mdx are shown. 
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K 


aVo 


a 


a 2 cj 


Srot 


a/r 


a/n 


ro/n 


0.156 


0.6371(20) 


0.2911(19) 


5.748 xl0~ 2 (45) 


0.3349(64) 


0.2057(7) 


0.2848(9) 


1.3845(9) 


0.1565 


0.6407(19) 


0.2929(20) 


5.363 xl0~ 2 (43) 


0.3292(59) 


0.1988(7) 


0.2754(8) 


1.3854(9) 


0.15675 


0.6457(17) 


0.2958(16) 


5.035xl0- 2 (34) 


0.3253(61) 


0.1928(6) 


0.2674(8) 


1.3867(7) 


0.157 


0.6463(18) 


0.2955(18) 


4.892xl0~ 2 (38) 


0.3218(64) 


0.1900(6) 


0.2635(8) 


1.3866(8) 


0.15725 


0.6507(16) 


0.3010(17) 


4.577 xl0~ 2 (34) 


0.3291(58) 


0.1842(6) 


0.2559(8) 


1.3892(8) 


0.1575 


0.6524(17) 


0.3007(16) 


4.416xl0^ 2 (36) 


0.3221(56) 


0.1809(7) 


0.2513(8) 


1.3890(8) 


0.15775 


0.6559(16) 


0.3022(16) 


4.186xl0" 2 (37) 


0.3184(55) 


0.1762(7) 


0.2449(9) 


1.3898(8) 


0.158 


0.6560(15) 


0.3013(16) 


3.936xl0~ 2 (34) 


0.3173(53) 


0.1708(7) 


0.2374(9) 


1.3894(8) 



TABLE I: The four fit parameters of the static potential and the extracted values of afro, ajr\ and ro/n for different K". 
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FIG. 6: The parameter a versus two determinations of am^ (corresponding to whether the lattice pion mass used for the quark mass deter- 
mination was taken from the PP or the AA correlator). The figure also shows the straight line fit that incorporates the lowest four am q values. 
The notation PP, AA and AP is explained in Ref. [ 10 1 and briefly in Sec. [TT] 



of the data analysis and as a result our conclusions do not depend on the particular values used. 
The values of the fit parameters aVo, a, a = a 2 a and 5rot are presented in TablepJ 



V. SOMMER SCALE 



After correcting for finite lattice effects, we compute the derivative of the corrected (continuum) potential aV c given by the 
first three terms of Eq. ( |4.2| > 

R 2 4-aV c = aR 2 -a . (5.1) 
dR 




FIG. 7: The parameter a<7 1//2 vesus am q . The figure also shows the fits for the lowest five am q values as explained in the text. 
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FIG. 8: The fit parameters oVq and 5rot versus two different determinations of am* p 



Phenomenological studies of the static potentials yield [8] r 2 *E | r=ro = 1 .65 at ro : 
as the Sommer scale. In general, there can be many such scales r c (in fm) defined by 



= ^ 



: 0.49 fm. The distance scale ro is known 



(5.2) 



where jV c is a dimensionless number like 1.65 and r c is a corresponding distance scale in fm obtained from the static potential. 
In our case, using Eq. ( |5.1| l, this means 



a 

I'r 



1 



(5.3) 



Thus if we can determine the potential accurately in the intermediate range around r c fm, we can determine the lattice scale. 

Since we have used R m [ n = v2 ~ 1.4 and R msx = 3V5 ^6.7 as the fit range for the static potential, the Sommer parameter 
ro ~ 0.49 fm (corresponding to Rq = ro/a ~ 6) is barely within this range, we have also used another Sommer scale r\ such that 
jV\ = 1.0 [7| (corresponding to R\ = r\/a~ 4.4). However, we have not observed any noticeable improvement in the accuracy 
of the results obtained using the scale r\ , because our fits describe the corrected potential data very accurately much beyond 
^max as already observed in Sec. IV 



For any other potential-derived distance scales like r\ (in fm) which are not known phenomenologically as accurately as ro, 
the strategy is to compute the ratio ro jr\ accurately on the lattice and determine r\ from the ratio. 



VI. DEPENDENCE ON am q 



Figs.|]|7]and[8]show respectively a, a 1 / 2 = aa 1 / 2 , aVo and 5rot versus quark mass (derived from PCAC) in lattice units at 
all the eight values of K. These figures show that while a and aVo show only about 3% change for the entire range of am q from 
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0.07 to below 0.015, ac 1 / 2 go through a change of about 20% in the same range. Table[l]which has a — a 2 a also shows about 
30% change in this range. The values of the coefficient of the 1 /^-correction term, 5rot, shown in the lower part of Fig. [8] are 
quite close to the a values and are also similarly weakly-dependent on am q . 

Moreover, Fig. [6]shows that, for small enough am q (am q < 0.035 with our data), a is independent of am q . That a is weakly 
dependent on am q and is roughly constant within our accuracy is also seen in Fig. [9] which shows that act 1 / 2 is proportional to 



a/ro and a/r\ at least for small enough am q and from Eq. (5.3 1 the proportionality constant is 



Kc = V^-a (6-1) 

with JV C = 1.65 and 1.0 respectively for a/ro and a/r\. From Fig. [9] we compute a = 0.3012(16) (from the a/ro dependence) 
and a = 0.3013(8) (from the ajr\ dependence) for am q < 0.035. This coincides with the value obtained from Fig. [6] viz., 
a = 0.3013(8) again for am q < 0.035. 

Fig. [7]shows ad 1 / 2 to be linearly dependent on am q with a positive intercept on the ad 1 / 2 axis: 

aa l/2 = C l +C 2 am q (6.2) 

The fits use the lightest 5 quark masses corresponding to am q < 0.04 and C\ and C2 are dimensionless constants. 

In the past, CP-PACS Collaboration tried non-linear chiral extrapolation for a 2 a in Ref. l24l and linear extrapolation of 
ad 1 / 2 later in Ref. [25|. We note that, UKQCD (26), CP-PACS El and SESAM and T^L collaborations JH have found 
am^ -independence of a in the small am q region for Wilson or 0(a) improved Wilson fermions. 



The dependence of a/ro and a/r\ on am q is shown in Fig. 10 Consistent with the am q -independence of a and the linear 



dependence of ac 1 / 2 , both a/ro and a/r\ depend linearly on am q at least for small enough am q . The linear Fits 



a 



= A C + B C amq (6.3) 
r c 

are done for the lightest 4 quark masses (corresponding to am q < 0.035). The dimensionless constants A c and B c are consistent 
with the a and the ao l l 2 fits, i.e., A c = C\/K c and B c = C^/Kc within our accuracy. 



VII. INTERPRETATION OF THE am q DEPENDENCE OF a / r c 

In the previous section, we have presented our numerical evidence of am q -dependence of a/r c . For small am q < 0.035, we 
have shown that this dependence is linear and can be attributed to a similar linear dependence of ao 1 / 2 on am q while a does not 
appear to depend on am q within this range. 

Given that our simulations are done with unimproved Wilson gauge and fermion action, one may attribute all the observed 
am q -dependence to scaling violations [29]. However, all other serious investigations of lattice QCD with dynamical quarks 
with a variety of improved gauge and fermion actions have also observed a significant dependence of a/r c on am q . For an 
early summary of this dependence with improved and unimproved actions see Ref. 1 30 1 . SESAM and T^L collaborations 
1 28 1 employed naive Wilson gauge and fermion actions and UKQCD collaboration IF3T1 employed naive Wilson gauge and 
O(a) improved Wilson fermion actions whereas CP-PACS collaboration [25| used both improved gauge and Wilson fermion 
actions. JLQCD Collaboration |32j using Wilson gauge action and O(a) improved Wilson fermion action has also observed this 
effect. Same phenomenon was noticed in simulations with standard and improved staggered fermions Q [33] [34] [35) . Recent 
simulations employing domain wall fermions [ 12|, twisted mass fermions [36] and overlap fermions ll37l have also observed the 
phenomenon. 

Obviously then, one cannot consider this purely as a cutoff effect. Some part of the quark mass dependence of a/r c , cor- 
responding to relatively large am q is to be considered scaling violations, but for small enough am q the dependence should be 
accepted as a physical effect lf38l [39l . 

Once one acknowledges that, for small enough am q , the linear am q dependence of a/r c is a physical effect and not an artifact 
of the cutoff, the natural choice consistent with a mass-independent scheme is to assume that the scale a is constant but r c 
changes with quark mass. To determine the scale a, one then needs to make a chiral extrapolation of a/r c to the physical point 
where estimates on r c are available from experimental data on the interquark potential in heavy-onium systems. In this scenario, 



the scale a drops out of Eq. (6.3 1 because it is the same for all values of the variables \/r c and m q 
1 



£/ c + B c m q with A c — as/ c (s^ c : a constant with dimensions of mass) (7.1) 

r r 



showing 1 / r c to have a chiral behavior similar to m p . 
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0.25 




FIG. 9: The parameter acr 1 / 2 



versus a/ro and ajr\. The figures also show the fits with the lowest four data points as explained in the text. 
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FIG. 10: Ratios a/r c versus am q . The fits are done with the lowest five am q values. 



In this scenario, a 1 / 2 also has a similar behavior (from Eq. (6.2 1): 



(J 1 / 2 = +C2 m q with C\ = dial (^i : a constant with dimensions of mass) (7.2) 

In our recent paper [ 10] where we have investigated the pion and the rho masses and their decay constants as functions of the 
quark mass, we have assumed the above mass-independent scheme and the scale a is independent of quark masses. 

In fact, application of chiral perturbation theory (%PT) is untenable if the scale a is not taken independent of quark mass. The 
left hand side of a %PT equation always contains a physical quantity (like m 2 ) which is scale-independent. On the right hand 
side, there are low-energy constants and quark masses which are scale dependent. Obviously, a chiral extrapolation using such 
equations are only possible if all the quark masses are determined at the same scale a. 

We again stress that for the above interpretation of a mass-independent scale to work, am q has to be small enough. In our 
simulation with 2 degenerate flavors of sea quarks and with unimproved Wilson gauge and fermion actions we find that for 
am„ < 0.035, a is independent of am q and aff 1 / 2 is linear in am„, observations that lead to the linear am q -dependence of a/r c 
(Eq. (6.3 1). If, in addition, the scale a is taken as independent of m q , Eq. (6.3 1 naturally evolves into Eq. ( |7.1| which shows 
dependence of 1 / r c on quark mass m q as a physical effect. 

One can consider the other extreme, i.e., r c independent of quark mass and the scale a changing with quark mass. This is 
considered, for example, in Refs. [31, 40 1, [41 1 and [42|. In this approach, the am q dependence shows up only as a scaling 
violation in the form B3ll44l 

In = ^ a k(P) k +N am q + higher order scale violations, (7.3) 
where a^'s and N are numbers. The j3-dependence is a reflection of asymptotic scaling |9|. However, it needs to be pointed 



out that, although written above for In (r c /a), Eq. (7.3 1 can be written down for any ln(l /aji) where ji is a hadronic scale like 
m n etc. Moreover, Eq. (7.3 1 does not include any non-perturbative dependence of ji (e.g. quark mass dependence of hadronic 
masses). 

However, in this approach, it follows from Eq. p.3\ that for small enough am q the dimensionfull a 1 1 2 is independent of the 
quark mass (because a is independent of am q for small am q ), something that apriori looks implausible because the string tension 
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Extrapolation 


Chiral limit of a/ro 


Chiral limit of ajr\ 


am n from PP 


am n from AA 


am n from PP 


am n from AA 




am q p 


i 


< 


amf^ 


amj: p 




am^ p 


a/r c = A C +B c am q 


0.1616(13) 


0.1627(10) 


0.1620(13) 


0.1618(12) 


0.2246(18) 


0.2262(13) 


0.2253(17) 


0.2249(15) 


K c = >JJfc - a 

1/9 _ _ 

aa 1 =C\ +C2am q 
a l r c = % + % am q 


0.1605(12) 


0.1617(10) 


0.1609(12) 


0.1607(11) 


0.2230(17) 


0.2246(13) 


0.2235(17) 


0.2233(14) 


a/r c =P c + Q c (am lt ) 1 


0.1631(16) 


0.1632(16) 


0.2274(20) 


0.2274(20) 


a/r c = D c + E c (m n /m p )' 1 




0.1591(37) 




0.2213(49) 



TABLE II: The values of a/r c extracted in the chiral limit using four different extrapolations. No entries for the extrapolation with (m K /m p ) 2 
when m n is determined from the PP correlator (because of finite size effects on the PP correlator at the largest two K* values, as expalined in 
text). 



is nothing but the energy density of the field flux between the heavy quark-antiquark pair and is likely to depend on the sea quark 
masses. There are of course other problems associated with such a mass-dependent scheme, e.g., usage of %pT and matching 
lattice scheme to mass-independent schemes like MS. 

Our observation of am q -independence of a for am q < 0.035 indicates that scaling violations in our data are negligible for 
small am q ; however, Eq. (7.3 i suggests that scaling violations are always present even for small am q . As a result we do not find 
support of the second scenario in our data and consequently pursue the first scenario where the scale a is taken as a constant for 
a given j3 and is determined from a chiral extrapolation of our a/r c data to the physical point. 



VIII. CHIRAL EXTRAPOLATION 



On the lattice am q is the best regulator of chiral symmetry breaking and as such it is the best parameter to use for chiral 
extrapolation of a/r c . However, it is not a good parameter for an approach to the physical point because apriori one does not 
know the value of am q at the physical point although it is very close to the chiral limit. 

For the chiral extrapolation of a/r c to the physical point, we have used (am n ) 2 . Fig. [TT| plots both a/ro and ajr\ versus 
(am n ) 2 . Firstly, we emphasize that we prefer (am K ) 2 rather than (r c m K ) 2 , because r c has its own quark mass dependence. 
Secondly, for the fit we stick to the linear part of the dependence corresponding to small {am K ) 2 . No higher powers of (am K ) 2 
is entertained to fit all data because as much as possible we want to stay away from data points which may have some scaling 
violations. Unfortunately, for a quantity like (am K ) 2 there is significant finite size effect at the smallest values with our lattice 
volumes. According to the findings of Ref. [ 10] we dropped the lowest two values when fitting with (am n ) 2 determined from 
PP correlators, and dropped only the lowest point when fitting with (am K ) 2 determined from AA correlators. 

We obtain the scale a by solving the quadratic equation in a: 

^ = P C + Q c (amf) 2 (8.1) 

where P c and Q c are constants, and rf h and the values at the physical point. 

We acknowledge that there is some unavoidable mistake made by using (am K ) 2 (as opposed to using am q ) for the chiral 
extrapolation because in the chiral region the relation between (am K ) 2 and am q is linear only in the so-called leading order. We 
feel that using {am % ) 2 is still better than using (r c m K ) 2 for reasons stated above and definitely better than using (m n /m p ) 2 which 
we also use for a rough estimate of the chiral extrapolation in Fig. 12 The squared ratio (m K /m p ) 2 is generally taken as an 



estimate of the quark mass, but it is not linear in quark mass for any appreciable range of quark mass. For small enough quark 



masses, Fig. 12 shows approximate linear behavior for the smallest masses. We have done the scale determination using chiral 
extrapolation of a/r c with respect to (m^/nip) 2 with pion masses determined only from the AA correlator because only in that 
case we have three data points (excluding the lightest masses at K = 0.158 but including the data at K = 0.15775) for a straight 
line fit. As pointed out by our earlier work II 101 . there are significantly more finite size effects on the pion mass determined from 
the PP correlator and we have to drop the lightest two masses (corresponding to K = 0.158 and 0.15775) and as a consequence 
would be left with only two points for a linear chiral extrapolation in this case. Hence we do the chiral extrapolation of a/r c in 
dependence of (m K /m p ) 2 when the pion mass is determined only from the AA correlator. 

In Table |n] we show the a/r c values at the chiral limit obtained with extrapolations done using am q (first data row), (am K ) 2 
(third data row) and (m n /m p ) 2 (fourth data row). The second data row contains the chiral limit values of a/r c obtained from 




FIG. 12: Ratios a/r c versus (m^/mp) 2 . The fits are done with the three lowest (m n /mp) 2 free of finite size effect. 
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FIG. 13: The ratio ro/n versus {am n ) . The fit is done with the four lowest (am n ) values. 





(fl/ n )/(fl/r ) 








amjt from PP 


am jt from AA 


versus 


r /n = Kq/K x 






am^ p 


arn^ 


aniq P 


{am K ) 2 






1.3900(16) 


1.3899(14) 


1.3901(15) 


1.3899(14) 


1.3894(7) 


1.3894(4) 



TABLE III: The ratio ro/rj in the chiral limit extracted using different methods. 



the individual limits of a and ad 1 / 2 . The first and the second row values are consistent with each other showing that our 
inference of a being almost independent of am q and flff 1 / 2 linear in am q is correct. Comparison of the values in the first and the 
third data rows shows that the central values are about 1% off and they are consistent with each other within statistical errors. 
These consistency checks give credibility to the chiral extrapolation of a/r c with respect to (am n ) 2 . Only the fourth data row 
containing extrapolated values using (m n /m p ) 2 shows a deviation of about 3% from the values in the first data row and also 
exhibit significantly larger statistical errors. The first data row is in bold font to emphasize that the data entries in this row have 
the most reliable chiral limits. 



Extrapolation 
to the physical point 


a/ro fit 


a/r\ fit 


am n from PP 


am n from AA 


am n from PP 


am K from AA 


a (fm) 


a -1 (GeV) 


a (fin) 


a -1 (GeV) 


a (fm) 


a' 1 (GeV) 


a (fm) 


a -1 (GeV) 


a/r c = P c + Q c {am K ) 2 


0.08027(77) 


2.458(23) 


0.08032(76) 


2.457(23) 


0.08053(70) 


2.450(21) 


0.08053(71) 


2.450(22) 


a/r c =D c + E c (m lt /m p ) 2 






0.07865(170) 


2.509(54) 






0.07873(164) 


2.506(52) 



TABLE IV: The lattice scale a (fm) and a 1 (GeV) obtained with two different extrapolations of a/r c to the physical point. 
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Extrapolation to 
the physical point 


am n from PP 


am n from AA 


a (fm) 


a' 1 (GeV) 


a (fm) 


a' 1 (GeV) 


amp = F\ + Fiiam^) 1 


0.07932(135) 


2.488(41) 


0.07995(195) 


2.468(60) 



TABLE V: The lattice scale a (fm) and a (GeV) obtained with extrapolation of amp to the physical point. 
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FIG. 14: amp vs (am n ) 



The fits are done with the lowest five pion masses free of finite size effects 



Fig. 13 plotted against (am K ) shows that for (am K ) < 0.1 or am q < 0.035 the ratio rg/fi is independent of these quantities. 



The fitted constant value in the figure is 1.3894(7) (shown in bold in Table III to indicate that this is the value actually used 
to determine rj). This value is absolutely consistent with the ratio of the chiral limits a/r\ to a/ro (extrapolated with different 
quark mass evaluations) and also with the ratio Kq/K\ = ((1.65 — a)/(1.00— a)) 1 / 2 (with the value of the fitted a put in from 
Fig. |6j, as shown in Table III 

We have used ro = 0.49 fm in our analysis. We are aware that a few lattice groups have calculated ro from the low level 
splittings of the bottomonium system and those values are a few percent lower than the standard value used in this paper. In 
any case we have done our simulation with 2 degenerate sea quarks only and since we have taken the view in this paper that 
r c changes with quark mass, it is conceivable that r c may change somewhat when the number of flavor is changed. We take 
the viewpoint that given that ro may have some uncertainties, we try to minimize all other uncertainties regarding the scale 
determination. 

Given the ro/n value at the physical point, with ro = 0.49 fm, we get r\ = 0.3527(2) fm. We can now plug in the values of 
ro and r\ respectively in the values of a/ro and a/r\ at the physical point obtained from the various extrapolations of a/r c and 



get the scale a. Table IV lists the values of a in fm and a in GeV obtained with the two methods employed. We notice that the 
scales obtained from (am H: ) 2 -extrapolation has very accurate values with less than 1% errors (emphasized by bold fonts in the 
table) while those obtained from (m K /m p ) 2 has about 2% errors. Within error bars the values are consistent. 

In our earlier paper ifTUl we computed the rho meson mass in lattice units and it was shown to have a linear behavior in am q 
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with a positive intercept at am q = 0. In Fig. 14 we plot am p as a function of (am K ) 2 . The data is well fit by a linear ansatz: 

am p = F\ +F2(am K ) 2 . (8.2) 

The smallest two and the largest pion masses are excluded from the fit for the case of the pion mass determined from the PP 
correlator while for the fit with pion mass determined from the AA correlator we have excluded the largest two and the lightest 
pion mass. These exclusions are due to possible finite size effects on the smallest pion masses and possibility of being outside 
the chiral regime for the larger pion masses. The value of the intercept at (am K ) 2 — is consistent with the intercept at am q = 
in Ref. iflOll . Again, the fit can be looked upon as a quadratic equation in the scale a (in fm) at the physical point while the pion 
mass m^ 1 and the rho mass mP p b are entered in fm -1 . The scales obtained by solving the equation are independent of the static 
potential and scales determined therefrom and are listed in Table |V| Although with relatively larger errors (~ 2 — 2.5%), these 



scales are very close to our very accurate evaluations using the a/r c extrapolations given in Table IV 



IX. ESTIMATE OF a 1 ! 1 



We can now make estimates of the parameter (7 1 / 2 in physical dimensions by chirally extrapolating act 1 / 2 with (am K ) 2 using 
the scale determined by our accurate determinations from a / r c extrapolations (first data row of Table IV i and the pion mass m^ h 
at the physical point: 

fl((J ph ) l/2 = Gi+G2 ^ m ph^ 2 (91) 

i included for the fits are similar 
ig the superscript indicating the 

value at the physical point) corresponding to two evaluations of pion masses from the PP and the AA correlator and two values 
of the scale a from a/ro and a/r\ extrapolations. All these four values are extremely close to each other. We present the average 
and quote the largest error of the four: 



The data along with the fits to determine the constants Gi and G2 are shown in Fig. [15] The points 
to the am p - (am n ) 2 fits discussed before. From Eq. (9.1 1 we obtain four values of?/ 2 (droppin 



rV2 : 



: 465.5 ± 1.4 MeV. 



(9.2) 



Similarly, using the scale determined independent of the static potential (from am p - (am n ) fits), a ' can be independently 
determined from the linear behavior of the ratio m p /c^/ 2 with am p (shown in Fig. 16 1. These determinations of c 1 / 2 have 



somewhat larger statistical errors than above. We present the average with the largest error: 

a 1 1 2 =460.9 ±9.3 MeV. 



(9.3) 



Although the statistical error in Eq. (9.2 1 is surprisingly very s mall, obviously it does not represent all the errors associated 
with the evaluation of a 1 / 2 as the second evaluation given in Eq. (9.3 1 shows a significant systematic deviation of the central 
values. 



X. IMPLICATIONS FROM WEAK COUPLING PERTURBATION THEORY 



Once we have determined the parameters of the static potential, it is interesting to compare the extracted parameters with those 
expected from weak coupling perturbation theory, wherever appropriate. For example, one would like to see how the extracted 
self energy quantitatively differs from that calculated in one loop (tadpole improved |46|) lattice perturbation theory. Moreover, 
from the non-perturbatively extracted average value of the plaquette, one can extract the strong coupling constant a v at the scale 
3.41 /a according to the Lepage-Mackenzie scheme 11461 . 

In weak coupling perturbation theory, the expectation value of the Wilson loop is given by 

LimitT^oc(-) -ln(W(R,T)) = V(R) (10.1) 

where 

V(R) = Vcoul + Vidf. (10.2) 

The static Coulomb potential 

Koui = -Cf ^ (10.3) 

K 

and V se if is the static source self energy. 
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FIG. 15: a^/a vs («i) 2 . The fits are done with the four lowest pion masses free of finite size effects. 



A. Static source self energy 



In the continuum, to lowest order, 



V S eif = C F 4na s 



d 3 q 1 



(2^) 3 q 2 



(10.4) 



On the lattice 



aVseif = C F 4na s -3 £ 



l3 q % L sin 2 a^/2 ' 



(10.5) 



Using /3 



5.6 and |f 



= 2.9987 for a 16 3 lattice, we get aV se if=0.3409. Incorporating tadpole improvement 



9 .j £ ,sin c "?'/ 2 

utilising the average value of the plaqutte (day), g 2 ^ g 2 = = g 2 /0.5744, we get flV se if=0.593. This value although calculated 
at the lowest order (without any quark loops) may be compared with the value we get from the numerical fit to the Wilson loop 
data, aVo ranges between 0.63 and 0.66 for the range for am q explored. 

The difference might be due to higher order corrections and/or nonperturbative contributions. 

We note that in perturbation theory, the strength of the static Coulomb potential and the static source self energy are both 
given by the strong coupling constant a s . In higher order of perturbation theory [451 the static source self energy can have 
am q -dependence, the self energy increasing with am q decreasing. As pointed out earlier and as shown in Fig. [8j aVo has a weak 
dependence on am q with the trend suggested by Ref. [45 1, but, however, at our smallest quark masses, it approximately saturates. 
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FIG. 16: The ratio m p / ^fa versus am p . The fit is done with the lowest five data points excluding the lowest point. 



B. Extraction of the strong coupling constant from the plaquette 

According to Lepage and Mackenzie l46ll . the strong coupling constant a v (q) at momentum scale q = 3.41 ja is defined via 



V(q) = -C F An , 
= -C F An ' 



0£o 

a v (q) 



/ (11 n 2 

1 + Oq . 3 — In 



An 



a 2 q 2 



And 



r 



where d = d g + d f = 0.37428 -2.x 0.00426 = 0.36576, C F = 4/3 and n/ is the number of flavors. 
Solving for the bare coupling 0£q = g 2 /An from 



a^) = ao(l + Oq C) with C = ^ . 3 J/ ln 



An 



a 2 q 2 



And 



we get 



ao = 0£v(l -C a v ) 

From perturbation theory [47 1, for the logarithm of the average Plaqutte (D av ) one gets 



In D av = cig 1 +[c 2 +-c\) g 4 



(10.6) 



(10.7) 



(10.8) 



(10.9) 



23 




FIG. 17: Average value of the plaquette versus 1/fC. 



with ci = 1 .3 and c 2 = cf + c{ = 0.03391 - .003696 = 0.0302. Thus 

-lnD av = ^47rao + 0.0858 x \6n 2 al 



= -Ana*, [l-Ov(C- 0.0858 x \2n)\. (10.10) 
From the measurement of the average plaquette (D av ) at each K and extrapolating to K c one can determine the strong coupling 



constant at a given scale. As shown in Fig. 17 omitting the data points at Jf=0.156 and 0.158 we perform a linear fit of the 
average plaquette and find □ av = 0.57438 at K c = 0.15857. Using this value, we find a v (3.41/a) = 0.167. If we run down the 
coupling to the scale I /a using two loop renormalization group formula, we get a v (l /a) — 0.257. In the same convention, the 
coefficient of the 1 /R potential from the fit (see Fig. |o]l a — \ x 0.30 = 0.225. 



XI. CONCLUSIONS 



Understanding the dependence of the static potential and the quantities derived from it, in particular a/r c , on the sea quark mass 
is mandatory to set the lattice scale using the potential and is also important for chiral extrapolation of hadronic observables. At 
present there is very little theoretical understanding of the sea quark mass dependence of the parameters a 1 / 2 and l/r c . Accurate 
numerical evidence at small am q of the observed dependence may be useful for theoretical understanding. 

Interestingly this dependence has been observed by all lattice QCD simulations even including the improved actions like 
overlap and domain wall fermions along with improved gauge actions. As a result it is difficult to consider this as a cut-off effect 
purely. 

Our approach has been to investigate the issue with an action with the most lattice artifacts, i.e., the standard Wilson fermion 
and gauge action, but with (i) a large enough j3 (=5.6) so that the scale a is small enough, and (ii) a large set (eight values) of 
the fermionic hopping parameter K for sea quarks corresponding to PCAC quark masses in lattice units am q from around 0.07 to 
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less than 0.015. This allows us to look for lattice artifacts or scale-violations quantitatively in terms of a reasonably large range 
of am q . 

What we have found numerically is that for small am q corresponding to am q < 0.035, the quantities related to the static 
potential have specific orderly behavior. With the usual Cornell potential parameterization, we find that for am q < 0.035, the 
parameter a (coefficient of the 1 /R term) is independent of am q while ao l l 2 (coefficient of the linear R term) depends on am q 
linearly with a positive intercept at am q = 0, resulting in a linear am q dependence of the quantities a/r c where r c is a Sommer- 
type scale parameter. We have taken utmost care in all aspects of the analysis to come to the above behaviors of quantities 
numerically, e.g., with regard to optimum smearing level for the gauge configurations at each K, correction for 1 /R on the finite 
lattice, fitting range [r m i n ,r max ] of Wilson loop data, fitting range [/?mi n , /?max] of the static potential data etc and we believe 
that these conclusions are independent of choice of parameters of the analysis at least qualitatively and to our precision even 
quantitatively. 

With the above, now if we accept a mass-independent scheme, the above linear am q dependence of acr 1 / 2 and a/r c naturally 
translates into a physical linear m q dependence of a 1 / 2 and l/r c . These dimensionful quantities are then very similar to m p as 
far as m q dependence is concerned. 

Once m q dependence is taken as a physical effect, in the mass independent scheme there is then conceptually no problem 
taking a chiral extrapolation of the numerical data of ad 1 / 2 and a/r c to the physical point for a dimensionful value of a 1 / 2 and 
the scale a. 

We have exercised care also in the chiral extrapolation. The quark mass in lattice unit, viz., am q , is the best quantity for a 
chiral extrapolation, but it is not suitable for an extrapolation to the physical point. We have first used extrapolation with respect 
to (am n ) 2 to make sure that we get the same limits at the chiral point, i.e., am q = or (am K ) 2 = 0. This was easily achieved 
once large (am^) 2 points and also (with knowledge of finite size effect from our previous work [10]) the smallest one or two 
(am n ) 2 point(s) were omitted from the fits. We stress that we stick to fits only with linear power of (am n ) 2 and we prefer (am K ) 2 
to (r c m n ) 2 or to (m^/mp) 2 , because (i) r c has its own quark mass dependence and (ii) (m % jmp) 2 is only approximately linear in 
am q even for small quark masses. 

We obtain an accurate determination of the scale by solving a quadratic equation in the scale a (in fm) resulting from the linear 
dependence of a/r c on (a»%) 2 and putting in the values of rf h and m^ 1 at the physical point. We quote with ~ 1% error: 

a = 0.08041 (12) (77) fm, a' 1 = 2.454(4) (23) GeV. (11.1) 

The first and the second errors are respectively the systematic and the statistical errors both of which are estimated conservatively. 



Out of the four jackknife statistical errors shown in the first data row of Table IV we quote in the above only the largest error. 
Also the systematic error is estimated by halving the systematic difference between the scale determinations from a/ro and ajr\. 

From (m K /m p ) 2 extrapolations also, we have determined the scale which is consistent with our accurate determination above 
but has large (> 2%) errors. 

In order to have an independent check on the scale, we have determined the scale also from a linear am p — (am n ) 2 extrapo- 
lation, a method which is fully independent of the static potential and the quantities derived from it. It is very satisfying to find 
the scale obtained this way comes within 1% of our accuarate determination achieved with the extrapolation of a/r c . Errors are, 
however, large (~ 2%) in this case. 

In our determination of the lattice scale from the static potential, we have assumed ro = 0.49 fm. There have been a few 
determinations of ro from low level energy splittings of the heavy-onium systems and these values are a few percent lower than 
the value we have used. There have also been suggestions in the literature about the uncertainty of the value of ro because it is 
not a quantity directly measurable from experiments. We like to mention that the hadronic masses also have some uncertainty 
in their values and like ro they are also likely to change with the number of flavors. Our approach in this paper has been that 
given the uncertainty in the value of ro, we wanted to reduce the uncertainty in the rest of the determination. In addition to 
the scale determination, the quark mass dependence of a/r c is also a physics issue that one needs to understand. Anyway, our 
scale determination independent of the potential shows that the scale obtained is extremely close to the value obtained from the 
potential. This indicates that the value of ro (= 0.49 fm) used in our analysis may not suffer from major uncertainties. However, 
we should keep in mind that our analysis is done with 2 flavors of sea quarks. 
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